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Abstract 

We seek to determine the mechanism of like-charge attraction by measuring the temperature 
dependence of critical divalent counterion concentration (Cc) for the aggregation of fd viruses. We 
find that an increase in temperature causes Cc to decrease, primarily due to a decrease in the 
dielectric constant (e) of the solvent. At a constant e, Cc is found to increase as the temperature 
increases. The effects of T and e on Cc can be combined to that of one parameter: Bjerrum length 
{Ib)- Cc decreases exponentially as Ib increases, suggesting that entropic effect of counterions plays 
an important role at the onset of bundle formation. 

PACS numbers: 87.15.Nii, 87.15.Tt, 87.16.Ac,87.16.Ka 
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Multivalent counterion induced aggregation of polyelectrol^es^ has been observed in a 

During the past 




variety of systems such as DNA, F-actin and bacteriophages [1|, 

decades, there have been intensive theoretical investigations to elucjxlate^he mechanism of 
counterion induced attractive interaction between polyelectrolytes 

n 

Il3l |. Although other types of interactions such as hydration 
contributing factors, it has been estabhshed that the electrostatic between the polyelec- 
..ol^e. and the. co„e.ated cou„te.„„s B B 3 B H H H is the .ajo. ca„.e of 
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like-charge attraction. 

There are, however, two possible ways, in which the condensed counterions on different 
polyelectrolytes correlate with each other. Thermal fluctuations create transient regions of 
high or low counterion density along the polyelectrolytes, which are typically simplified to 
be charged lines each with a thin layer of condensed counterions When two parallel 
polyelectrolytes get close to each other, an attractive force is induced by their transient 
complementary counterion density profiles. The attractive force induced by the long wave- 
length counterion fluctuations has been predicted to be proportional to temperature using 
the mean field theory Alternatively, the counterions may correlate with each other in 
their positions on the surfaces of polyelectrolytes. A representative picture of the attrac- 
tive interaction induced by positional correlations of counterions is provided by the Wigner 
crystal model 0. In this model, condensed counterions form Wigner crystals on the poly- 
electrolyte's surface at T = K. Cross correlation of counterions occurs when the distance 
between two polyelectrolytes decreases to the lattice constant of the Wigner crystals. It is 
the cohesive energy of the Wigner crystals that results in attractive interaction and aggrega- 



tion of like-charged rods 
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. Since thermal fluctuation diminishes the structural order of 



counterions, this model predicts a stronger attractive interaction at a lower temperature . 
Recent computer simulation results suggest that the Wigner crystal model captures the 



physics of like-charge attraction better than the thermal fluctuation model p], since it takes 
into account the strong correlations between individual counterions. However, the Wigner 
crystal model is a zero temperature approximation. It is uncertain whether the model 
works under high temperatures, even though it has been suggested in the theoretical work 



that the short range order of counterions survives at high temperature [6|, 3, ISl • Since the 
two theoretical models predict opposite temperature effects on attractive interactions, the 
aggregation of polyelectrolytes under different temperatures is expected to directly reveal 
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the dominant mechanism of hke-charge attraction. 

In this paper, we present the results of temperature effects on critical concentration of 
counterions Cc, using MgClg and CaCl2 to induce bundle formation of bacteriophage fd. We 
found that Cc increases with increasing temperature under a certain dielectric constant. The 
effect of the dielectric constant on Cc has been measured as well. Under a fixed temperature, 
Cc decreases with decreasing e. Noting the similar effects of T and e on threshold concen- 
tration, we speculate that Cc is actually determined by the Bjerrum length {Ib)- The result 
that Cc doesn't change with temperature at a constant Ib validates our speculation. The 
thermodynamics for aggregation of fd virus is investigated using the result that Cc drops 
exponentially as a function of Ib- The effect of hydration of counterions can be assessed 
from the prefactor in the exponential expression. Finally, based on the result that stronger 
counterion correlations lead to lower threshold concentration, we conclude that our results 
support the Wigner Crystal model. 

Bacteriophage fd is a rod-like polyelectrolyte of approximately 880 nm in length, and 
6.6 nm in diameter. There are approximately 2,700 copies of coat protein arranged crys- 
tallographically on the virus surface. At neutral pH, each coat protein contributes four 
net negative charges on the surface of a virus. Thus, a virus has a linear charge density 
of approximately 12.5 e/nm. When a virus is simplified to a linear array of charges, the 
charge spacing would be 6 = 0.8 A. The fd viruses are highly charged so that their latteral 
aggregation can be induced by divalent counterions such as Mg^^ and Ca^^ HQ. The 



viruses survive at temperatures up to 90 °C {20! |2l|. Therefore, within our experimental 
range (T < 50 °C), the results are not expected to be due to the changes in the physical 
parameters of virus. 

Aggregation of fd virus is detected by measuring the scattering light intensity at a fixed 
angle of 90°, with a PERKIN ELMER LS-5 luminescence spectrometer. 800 /il 0.1 mg/ml 
bacteriophage fd was added to a rectangular cuvette of 10 mm path length and 5 mm width. 
Scattering intensity was measured when the solution reached its steady state following the 
addition of a stock solution of concentrated CaCl2 or MgCl2. Sample temperature was con- 
trolled using a ISO TEMP 1006S water bath (Fisher Sci Inc.) connected to the sample 
holder. When increasing the counterion concentration, an abrupt increase in the light scat- 
tering intensity was noted at a critical concentration Cc. The Cc was defined at the divalent 
concentration where at least a 20-fold increase in the total scattering intensity was observed. 
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To find out the temperature effect on tliresliold concentration, we first measured tlie 
tfiresliold concentration under different temperatures without other manipulation of solution 
properties. Next, the effect of solution dielectric constant on Cc was measured under a 
constant temperature of 20 °C. Finally, Cc was measured at different temperatures with the 
dielectric constant fixed. 

Fig ]l(a)| shows threshold concentrations of MgCl2 and CaCl2 as functions of the sample 
temperature. Cc for both MgCl2 and CaCl2 decreases with increasing temperature. However, 
the temperature changes of a solution also results in changes of its dielectric constant. The 
dielectric constant of the solution, which is approximately equal to that of water, under 
different temperatures are calculated using e{T) = a + bT + cT^ + dT^, where a , b , c , d 
are empirical constants [22]. It is shown in Fig |l(a)| that the calculated e of water decreases 
from 83.9 to 73.2 when temperature increases from 10 to 40 °C. 

Changing the dielectric constant dramatically affects Cc. In Fig ]l(b)[ Cc is measured 
under different ethanol concentrations at T = 20 °C. Adding ethanol to a solution effectively 
changes its dielectric constant, e for a mixture of ethanol and water is determined by 
e(T) = (1 — Ce)ew(T) + CeEeiT) [2^, where ew(T) and ee(T) are dielectric constants of 
water and ethanol at temperature T respectively, and Cg is the volume fraction of ethanol 
in the solution. At 20 °C, = 80.37 and Ee = 25.00. According to the relation between e 
and Ce, the dielectric constant of a solution drops when its ethanol concentration increases 
(Fig |l(b)1 . The data in Fig jl(b)| show that a decrease in e leads to decrease of Cc. Similar 



effect of ethanol on critica 
condensation experiments 



;rivalent counterion concentration has been observed in DNA 
^. The relationship between Cc and e indicates that the like- 
charge attraction originates from electrostatic correlations between counterions. Since the 
strength of counterion correlations is proportional to 1/e, a smaller e leads to a stronger 
counterion correlation, indicated by a lower Cc. 

Since Cc decreases upon the decrease of e, the dependence of Cc as a function of temper- 
ature under a fixed dielectric constant is required to manifest the pure temperature effect. 
In FigEl Cc increases as T increases under fixed dielectric constant, e of the sample is 
tuned to be constant under different temperatures by the addition of appropriate amounts 
of ethanol. At e = 73.2, Cc for MgCl2 increases from 10 mM to 70 mM as we raise the 
temperature from 10 °C to 40 °C. Hence, the drop of Cc in Fig |l(a)| is actually the net effect 
of an increase in temperature and a decrease of the dielectric constant, which is induced by 
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the rising temperature. The curve for e = 69.9 is lower than that of e = 73.3. This is to be 
expected, because lower dielectric constants lead to lower threshold concentrations. 

A decrease in either e or temperature leads to a lower threshold concentration. We 
speculate that the effect of temperature and dielectric constant can be combined to that of 
one parameter: the Bjerrum length Ib = e"^ /AireeokBT, which is often applied to quantify the 
importance of counterion correlations through a coupling parameter S = 2iTz^l'^a, where z is 
the valence of counterions and a represents the surface charge density of polyelectrolytes 



25|. At strong correlation limit, i.e. S ^ 1, the condensed counterions form two dimensional 
(2-D) strongly-correlated liquid or even Wigner crystals at low temperatures. Cc is indeed 
a function of Ib, since it doesn't change with temperature at a constant Bjerrum length 
(See Fig 121). More importantly, a lower Cc is observed at a lower Bjerrum length. Thus, a 
stronger correlation between counterions leads to a smaller critical counterion concentration, 
since larger Ib values indicate stronger counterion correlation. 

The threshold concentration is quite sensitive to the change in Bjerrum length. In FiglHl 
Cc of MgCl2 and CaCl2 decrease exponentially as a function of Ib- The logarithm of Cc fits 
to a linear function of Bjerrum length: In Cc = —AIb + B , with A = 1.89 A ^ , B = 17.6 for 
MgCl2 and A = 2.09 A ^ , B = 18.5 for CaCl2, respectively. The linear relationship between 
InCc and Ib indicates that the counterion entropy plays an important role in the process 
of aggregation. The difference in entropy for the counterions in a bundle and those in the 
bulk medium is approximately InCc/Co, where Co is the counterion concentration in the 
condensation layer. This entropy loss of the condensed counterions causes an increase in the 
energy barrier which has to be overcome by the attractive interaction induced by electrostatic 
correlations. Since the energy gain from electrostatic interactions is proportional to Ib, a 
linear relation between In Cc and Bjerrum length is expected. 



Based on the mean field theory approximation, t 



le Oosawa theory predicts that the force 



per unit length between two polyelectrolytes as 

1 z^e 



f{R) = kBT{ 



z^RIb il + z^e)R^ 

where z is the valence of counterions, R is the distance between the long axes of two parallel 
polyelectrolytes, and ^ = Z^/fe is the Manning parameter. For highly charged polyelec- 
trolytes, ,^ ^ 1, f{R) is attractive when the distance R between two charged lines is less 
than z^Ib- To induce aggregation of fd virus with divalent counterions {z = 2) at a large 



Ib, say /fi = 9 A, the distance between two neighboring viruses must be smaller than 4 nm. 
However, the distance between the long axes of two parallel viruses could not be smaller than 
the diameter of a virus, which is approximately 6.6 nm. Therefore, it would be impossible for 
divalent counterions to induce aggregations of fd virus solely based on this treatment. The 
Oosawa theory fails because of its mean field nature and the simplification of polyelectrolytes 
as charged lines. 

The Wigner crystal model takes into account the strong correlation between counteri- 
ons. Before aggregation, each fd virus is surrounded by a layer of condensed counterions, 
which partially compensates the charge of the fd virus. As a result of the decrease in fd 
surface charge, the two viruses can move close to each other by Brownian motion, since the 
Coulomb repulsion is reduced. When the distance between two viruses is small enough, the 
counterions on different viruses start to correlate with each other. This correlation creates 
an attractive interaction between two viruses. Each counterion gains an energy of Gcorr from 
this attractive interaction. Since S ^ 14.8 for fd virus at room temperature and e = 80, 



the condensed counterions can be treated as a strongly correlated liquid, and Gcorr is t 



estimated as the energy of interaction between a counterion and its Wigner-Seitz cell 



then 

0. 



Gcorr = ^{zb/rf%kBT (1) 

where b is the charge spacing on the polyelectrolyte when considering it as a charged line, 
and r is the radius of a virus. Gcorr is estimated to be 0.13 IsksT per counterion within a 
bundle of fd viruses using 6 = 0.8 A and r = 33 A. 

The strong correlation also drives more counterions from the bulk solution into the space 
between two viruses ^,0]. These extra counterions roughly neutralize the residual charge 
of fd viruses, thus diminishing the Coulomb repulsion. However, the entropy loss of these 
counterions incurs an energy penalty. Assume that a fraction 7 of all counterions in the 
condensation layer are the extra ones, the average energy penalty per counterion is 

G« = 7fcBTlnCo/C (2) 

Aggregation of fd viruses occurs when the energy gain is larger than the energy penalty, 
Gcorr — Gr > 0. Using Eqs.(IT)) and the critical concentration Cc can be calculated as 

7/41^1/4 

lna = -^;^/B + lnCo (3) 
6 



From Eq.Q, 7 can be estimated by setting the fitting parameter A = ■ Based on 

the values of A = 1.89 A"^ and 2.09 A~\ we obtain 7 ^ 6% and 7% for CaCla and MgCla, 
respectively. When calculated from the Manning counterion condensation theory, approx- 
imately 94% of the surface charge of a virus is neutralized by divalent counterions at the 
onset of aggregation. These predictions suggest that the viruses are totally neutralized when 
they aggregate. 

The second term in Eq.Q contributes to the fitting parameter B. The counterion con- 
centration in the condensation layer has been measured to be a few molars [2^. Using 
Co = 3 molar for both Mg^+ and Ca^+, the parameter B is estimated to be 8.01, which 
is much smaller than the fitting result. The discrepancy likely originates in the change in 
counterion hydration state. A counterion loses its waters of hydration when it resides in the 
condensation layer. The ^^Mg NMR measurements of F-actin solutions indicate that only a 
small fraction of the condensed Mg^"*" lose waters of hydration 23|. If a small portion, say 
10%, of the extra counterions lose their waters of hydration, an energy penalty of approxi- 
mately lOkBT is expected for each extra counterion, since the energy penalty for a Mg^+ or 
Ca^+ ion to lose all of its waters of hydration is on the order of lOOfc^^ [28!|. When taking 
into account the hydration effect, Gr = ■-ykBT(\nCo/C + 10). Thus, S = InCo + 10 ~ 18, 
which is close to the fitting results. This estimate suggests that partial loss of hydration 
contributes to the change in Gibbs free energy, thereby affecting the onset concentration of 
multivalent counterions for bundle formation. 

In the DNA condensation experiments 12| , the threshold concentration of Mn^"*" ions to 



induce DNA condensation has also been reported to decrease with increasing temperature. 
The phenomenon is similar to the result in Fig|l(a)j When taking into account the variation 
in solution dielectric constant, the data in Ref. (l^ indicate an exponential decay of Cc with 
increasing Ib as well. Instead of discussing the effect of variations in dielectric constant, 
the data was interpreted by assuming that the attractive force originates in the release of 



ordered water from the polyelectrolyte surfaces [1^ It is interesting that the analysis 
based on the water exclusion also predicts an extra binding of counterions in the process of 
DNA condensation Q]. 

In summary, we have studied the thermodynamic effect of like-charge attraction by mea- 
suring the temperature effects on threshold concentrations for divalent counterions to induce 
bundle formation of fd viruses. Increases in both the temperature and the dielectric constant 



7 



of the solvent lead to lower critical counterion concentration. We found that Cc varies as 
a function of Bjerrum length, which combines the effects of temperature and the dielectric 
constant. The Oosawa model fails to predict the aggregations of fd virus because it simplifies 
the polyelectrolyte rods to charged lines. The linear relation between InCc and Ib is well 
interpreted using the Wigner crystal model, supporting the argument that it is applicable 
at finite temperatures. 

The authors appreciate the helpful discussions with Prof. B. I. Shklovskii, Dr. Thomas 
Angelini, and Dr. Don Rau. The work is supported by NSF DMR0405156 and NIH ROl 
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FIG. 1: Increase of temperature leads to decrease of dielectric constant (e) of the solvent, thus 
causes Cc to decrease, (a) When the e of the solvent is not fixed, Cc of both MgCl2 (solid 
triangles) and CaCl2 (solid circles) decreases with increasing temperature, e (solid squares) drops 
as temperature increases, (b) Cc of MgCl2 (circles) decreases with the amount of ethanol at 20 °C. 
Addition of ethanol leads to decrease of e (squares) . 
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FIG. 2: Effect of temperature on Cc under fixed dielectric constants and fixed Ib- At e = 73.2 
(circles) and e = 69.9 (squares), Cc increases with increasing temperature. When the solution 
dielectric constant is adjusted to fix Ib, Cc doesn't change with temperature . At fixed Bjerrum 
lengths, Cc is 45 mM with Ib = 7.28 (up triangles), and 40 mM with Ib = 7.38 (down triangles). 
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FIG. 3: Cc as a function oUb- Results for MgCl2 (filled triangles) and CaCl2 (filled circles) indicate 
that Cc decreases exponentially with increasing Is- Dashed lines are results of exponential fittings 
to the experimental results. 
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